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We analyze consequences of dynamical modulations of Josephson current by non-polar lattice 
mode in the Josephson junction barrier. In the high T c junctions, the effect of such modulations can 
be anomalously strong due to the proximity of the insulating barrier to the superconducting state. 
Accordingly, the interconversion of sound (as well as other non-polar phonons) and the Josephson 
plasma oscillations mediated by stationary Josephson currents, which may be present in the junction 
due to various reasons, becomes possible. We suggest that this effect can be employed for imaging 
of the stationary Josephson currents. Estimates of the effect are given. 



I. INTRODUCTION 



The nature of the superconductivity in the high T c materials continues to be a subject of hot debates. Following 
the suggestion very significant advances have been made in revealing the dominant d-wave symmetry of the 
superconducting order parameter (OP) in these materials A direct consequence of this symmetry is the possibility 
of the 7T- Josephson coupling M. The existence of the 7r-type Josephson junction (J J) and feasibility of its controlled 
implementation in circuits have been firmly established |^|,^ . The most startling consequence of the 7r-type coupling is 
the spontaneous Josephson current which generates the magnetic flux characterized by the half of the standard unit. 
It should be noted that a strong interest to the unconventional JJ made of the high-T c materials is inspired, besides 
purely academic attention, by the very recent proposal (^] of employing such J J as a qubit for quantum computations. 

Another intriguing consequence of the non-traditional symmetry of the OP is the time-reversal symmetry breaking 
(TRSB) at the grain boundaries Q. An automatic consequence of such states are spontaneous currents leading to 
fractional vortices, characterized by the magnetic flux which is neither the half nor the whole unit. It is worth noting 
that spontaneous fractional vortices have been detected in recent years [[tJ. Their exact nature, however, is still 
alluding the explanation. 

Detection of single vortices is usually achieved by means of the Scanning SQUID Microscopy technique (see in Rcf. 
H ) . It relies on measuring the total magnetic flux threading the pick up loop which scans the surface from the distance 
of few /j,m. Thus, only those spontaneous currents can be revealed which produce magnetic field normal to the surface. 
In fact, the resolution of the method || is limited by the pick up loop size (4-18/^m) as well as by the distance between 
the loop and the surface. Furthermore, the inductive interaction between the loop and the fluxes may distort the 
actual current distribution. In this regard we note that it would be desirable to have a technique capable of a direct 
high resolution imaging of the spatial distribution of the stationary (spontaneous) Josephson currents regardless of 
their orientation and of the scale of their spatial variations. 

In this paper, we discuss the effect which could be employed for the imaging of the Josephson currents with high 
resolution. This effect is based on the empirical observations of a very strong sensitivity of values of the Josephson 
critical current density J c in the grain boundary (GB) JJ to the misorientation angle between the grains ||. The 
nature of this effect remains very controversial. Many explanations have been proposed. One of them relies on the 
d-wave symmetry of the OP and faceting of the boundary, so that the macroscopic (average) critical current density 
is much smaller than the local critical current density |To| , which may change sign due to the 7r-type coupling. As 
suggested in Ref. an extremely strong suppression of J c is a result of mechanical deformations concentrated 
at the GB. It has been suggested that superconductivity is suppressed locally as long as the local mechanical strain 
becomes > 1% |ll]]. In accordance with the phenomenological approach jl^], the strain shifts the chemical potential 
of the JJ barrier, and this may strongly affect the Josephson coupling, as long as the barrier is characterized by the 
proximity to the metal-insulator transition. Many properties of the GB JJ have been successfully explained by this 
model. On the microscopic level, in Ref. (l3) it has been pointed out that small changes of the local structure, existing 
in the vicinity of the GB, should induce large variations of the valence of the copper ions, which control the free charge 
carriers concentration. 

Recently, the anomalous proximity effect has been induced by light in the Y-Ba-Cu-0 J J jlj]. This indicates that 
the JJ barrier can easily be perturbed (in this case, by light) into the superconducting state. As suggested in Ref. 
p5[ , this effect can naturally be explained if the electron system in the barrier is a sort of nearly superconducting 
quantum critical liquid, so that any external factor can easily shift its equilibrium. In fact, such a factor could be - 
light; mechanical deformations; optical displacements of atoms, etc. 
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In this paper, we pose a general question about mechanisms by which deformations (lattice displacements, in 
general) affect the Josephson coupling in the GB JJ made of the high-T c materials: 

Is the role of the deformations primary or secondary? 
If the effect of deformations is primary with respect to the electronic system, it can be revealed by means of applying 
reasonable deformations externally in the superconducting state and observing subsequent strong changes of J c . If it 
is secondary, the deformations may not affect directly the electronic system. These may predetermine the electronic 
properties during the formation of the GB JJ at high temperatures (way above the superconducting temperatures) 
through the processes of, e.g., oxygen depletion controlled by local mechanical stress and oxygen diffusion, so that no 
significant direct effect due to the external deformations should be expected to occur in the superconducting phase. 
In order to resolve this problem we note that, if the deformations play the primary role, a strong dynamical coupling 
should exist between sound and the Josephson phase. In contrast, if the deformations play a secondary role, no such 
a significant dynamical coupling should be observed. 

In this paper we will assume that the deformations play the primary role in accordance with the general concepts 
jn^plil ) and will concentrate on the dynamical implications of the strong sensitivity of J c to the mechanical and, 
possibly, optical (non-polar) deformations of the crystal lattice. Specifically, we will analyze a coupled dynamics of 
the phonons and the Josephson phase. We suggest that the coupling between the phase and the phonons should be 
anomalously large because of the proximity of the barrier to the superconducting state. In this phenomenological 
analysis, no quasi-particle effects are considered. This approximation can be well justified by the requirement that the 
phonon frequencies under consideration are much smaller than a typical superconducting gap. In this sense, we will be 
employing an adiabatic approximation in which the superconducting electronic ensemble adjusts itself momentarily 
to Josephson phase and the phonon variable(s). Thus, the following consideration is limited by the framework of the 
coupled dynamical Josephson equation |l6[ ] and equation of motion for the phonon variable. 

It is important to emphasize, that the first order coupling between small oscillations of the Josephson phase and 
sound (or optical non-polar phonons) can exist in a medium characterized by the inversion symmetry, if only some 
equilibrium Josephson currents are present. In other words, the linear coupling between the phase and the non-polar 
phonons can exist if and only if the time reversal symmetry is broken (either spontaneously or externally). This 
contrasts with the case of the polar phonons, which are coupled to the phase oscillations through the dynamical 
electric fields (due to the Josephson relation) regardless of the presence of the equilibrium currents. In this respect, 
we limit our consideration by the sound and the non-polar optical phonons, so that a mere presence of the linear 
coupling between the phase and the non-polar phonons is indicative of the stationary Josephson currents (we do not 
consider the case when no inversion symmetry exists). We believe that predictions following from our analysis provide 
an opportunity for experimental testing of the dynamical properties of the JJ barrier in the high-T c materials, and 
for answering fundamental questions about the nature of the Josephson effect in the high-T c materials. Furthermore, 
the discussed effect could be employed for the imaging of the stationary Josephson currents. This, in its turn, may 
shed light on the nature of the spontaneous Josephson currents and the TRSB in these materials. 

In Sec. II, we show how the linear coupling between the Josephson phase and the non-polar phonons is induced by 
the stationary Josephson currents. Then, the implications of this coupling are considered. Specifically, in Sec. Ill, the 
effect of the conversion of the non-polar phonons into the phase oscillations are discussed. In Sec. IV, we consider the 
induced infra-red (IR) activity of the non-polar optical phonons in the JJ. Then, in Sec. V, the effect of the coherent 
generation of sound in the JJ by the incoming microwave radiation is analyzed, and the possibility of employing this 
effect for the imaging of the currents is addressed. The role of external magnetic field in enhancing the interconversion 
is discussed in Sec. VI. Finally, the discussion and conclusion are given. 



II. COUPLING BETWEEN THE JOSEPHSON PHASE AND NON-POLAR PHONONS 

The Josephson current density J ((f) as a function of the phase difference tp can be sensitive to atomic displacements 
in the barrier. Let us assume that £(x, t) is some phonon variable. It can be the strain tensor Uij, or some other 
variable describing optical excitations. If this variable does not correspond to a polar (IR-active ) excitation, a direct 
coupling between £ and the phase oscillations is forbidden. In what follows, we will show how this non-polar phonon 
can become IR-active in the presence of the stationary Josephson currents, so that monitoring the IR absorption at 
the frequency of the crystal mode, which is known to be not IR-active in the absence of the static currents, would 
provide some information about these currents. 

In our approach, the Josephson equation for the phase ip(x,t) contains terms ~ £(x, t). If £, which modifies 
somehow the tunneling matrix element between the sides of the J J, is small, the Josephson coupling energy density 
(per unit area of the junction) Ej = Ej((p,£) can be expanded in £. The symmetry consideration may require 
that no linear term ~ £ is present in the ideal lattice. This is exactly the case for, e.g., the non-polar vector and 
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quadru-polar modes. If, however, £ is the strain, the scalar combination £ = Ujj (we employ the convention of 
the summation over the repeated indices) is allowed in the expansion. It should also be taken into account that the 
crystal structure close to the GB and at the GB can be highly distorted. For example, some stationary strain uf^ is 
inevitably present in any GB junction just to ensure matching of the sides. This, relaxes the symmetry requirement, 
so that the linear term becomes allowed. Thus, the expansion has a form 

Ej = Ef\^[l + b^,t) + o(e)}, (1) 

where Ej*\ip) denotes the Josephson energy density as a function of the phase in the absence of £; the coefficient 
b is taken as some constant. In fact, the observations of the very strong sensitivity of the critical current to the 

misorientation angle, which could be due to the local deformations arising in order to insure matching of the sides, 
justify the values as large as |6| « 10 2 . Indeed, it has been suggested that the deformations uq as small as 1% may 
cause the reduction of the current by almost an order of magnitude pT| . Thus, we take \b\ w l/u ~ 10 2 . As we 
mentioned above, the optical deformations of the lattice may produce a similar effect. In this case, the estimate for 
b can be obtained from the following considerations: optical deformation £ of 1% of the lattice cell (of the size 
a w 0.4nm) may produce the same effect as the strain, thus |6| « 10 2 /a. It should also be noted that the form (|l|) 
implies that the phonons do not modify the phase dependence of the JJ energy. 
The Josephson current, then, follows from eq.(|l]) as |l6[ ] 

J{<P,Q = ¥ ^ = M<P)[1 + &£(x,t) +o(e% (2) 

where Jo(<p) = {2e/Ti)dE^ / dtp is the Josephson current in the absence of £, with e > being the unit charge. 

The total Josephson energy of the system is obtained as the integral Hj — J cPxEj over the area of the JJ. The 
energy H% of the field £ can be represented in the spirit of a general Landau expansion with respect to the smallness 
of £ as 

H, = J d*x[^ + ^ + + o(£ 4 )], (3) 

where p, Uq, v g are phenomenological constants, so that p carries the meaning of some effective mass density. 
Accordingly, ujq and v g determine the dispersion law of small oscillations of £. 

The total energy H is the sum of H^+Hj and of the electro-magnetic energy Hem = J d x[CV /2+K,(V(p) /2] of 
the J J fl6fl , where C stands for the electric capacitance of the JJ per unit area; K — h 2 c 2 /(Wnde 2 ), with c, d standing 
for the speed of light, and the magnetic thickness of the JJ (given as the sum of the geometrical thickness db of 
the barrier and of the twice of the London penetration length Jig]), respectively. We take into account the Josephson 
relation 

V=^ (4) 
between ip(x, t) and the voltage V = V(x, t) across the JJ, and find 



H = Ht + Hem + Hj = J d'x[^- + ^ + '-fi^O 

(5) 

d 2 4j (^) 2 ^ 2 + + + m 

where the integrations J d 3 x and J d 2 x... are performed over the bulk and over the JJ barrier area, respectively. 
The above form yields equations of motion of the joint dynamics of (p and £. In this paper we will concentrate on 
the effect of small oscillations. 

The Josephson equation following from eq.(|^) is 

(^ + 7^-c 2 VV + *a/oM(1 + 6£) -0, f=||, (6) 



where all the variables are taken at the JJ plane; c = c/ \j4mdC stands for the effective speed of light in the J J; the 
term 799 describes the effect of dissipation in the resistively shunted J J |l6| ], and 7 = l/(p s C), with pj 1 being the 
shunt conductance per unit area of the J J. 
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The equation for £ also follows from eq.(||) as 

£ + ^ v 2 g VH + -Ef (<p)8(z) = 0, (7) 

where we have specified the JJ barrier as the plane z = 0. Note that, while the Laplace operator V = + in 

eq.@ acts along the J J plane, the Laplace operator V 2 = + + in eq.(|^) acts in the bulk. 

Equilibrium currents correspond to a non-trivial static solution of eqs.(^), (Q) which in the static case become 

- c 2 VVo + ^o(<A>)(l + = 0)) = 0, (8) 

wfco - «|V 2 £ + h -Ef{^)5(z) = 0, (9) 

where £o(z — 0) stands for the solution of eq.@ taken at the J J plane z = 0. In the long wave limit, the solution 
of eq.(^|) can easily be found and substituted into eq.(^), so that it acquires the form 

-c 2 V 2 ^ + ^Jo(^o)[l-54 0) (^o)] = 0, g= b (10) 

2p\0JQVg\ 

It can now be seen that the interaction with phonons modifies the equilibrium Josephson current, so that the effective 
current becomes J(y>) = Jo(</?)[l — gE^p{ip)\. The resulting current may significantly deviate from Jo(<p), if the 
coefficient g is large enough. This may occur if b is anomalously large, or the phonon mode exhibits soft-mode 
behavior. In this paper we will not consider such possibilities, and therefore it should be clearly shown that the 
dynamical effect under consideration does not imply strong deviations of the current-phase relation (CPR) from the 
bare one Jo (<£>), even though the values of b are taken large as discussed above (see below eq.(|l|)). To this end, we 
assume, for a while, that the current obeys the standard CPR Jo = J c smcp |l(|. Then, the effective current becomes 
J = J c (l — gfiJc/ (2e)) sin ip + (g?i J 2 /(4e)) sin 2ip. We require that the second term, which is ~ sin 2tp, is much smaller 
than the main one ~ sin ip, that is, 

_ _ gh.J c = b 2 hj c ^ 1 
4e 8ep\u!oVg\ 

Later we will show that for typical values of the parameters, g w 10~ 2 — 10 -3 , so that eq.(||) as well as the term 
~ £o m eq.(||) can be safely ignored. 

Here we do not specify mechanisms which may result in deviations of Jo(<p) from the standard CPR, and, as a 
possible consequence, in spontaneous Josephson currents (see, e.g., [^,0). Our goal here is studying implications of 
the presence of such currents on the process under consideration. In fact, the equilibrium currents may occur due to 
some trapped integer vortices or external magnetic field as well. 

Below we will consider coupling between small oscillations of (p, which constitute small deviations from the equi- 
librium, and £. We introduce small deviations ipfe, t) °f f from (po as 

tp{x,t)=<po{x)+i/>(x,t). (12) 

The variable £ should also be represented as the sum of the static solution £o of eq.(^|) and the oscillatory part. 
We will, however, ignore the static part, provided eq.(|i"l|) holds. 

The form (|l|) can be employed in eq.(|^) in order to obtain a linearized equation for ip. Before, however, we will 
show explicitly that, as long as the spontaneous currents are present, the variables ip and £ arc coupled linearly. 
Indeed, substituting (|l2|) into eq. (j^) and selecting the lowest order terms describing the interaction between ip and 
£, we find the interaction energy as 



/Tib 
d 2 x—J Q (tp )ip£, , 



(13) 



where £ is taken at the J J plane, and the definition of Jo(<p) (see below eq.(f|)) has been employed. Eq.(|l5|) indicates 
that the spontaneous (equilibrium) currents Jo((fo) induce the first order coupling between small oscillations ip and 
t 

It is worth noting that, the form ( |l3| ) of the coupling between the Josephson phase oscillations tp and a phonon 
variable £ is the lowest possible coupling as long as £ describes a non-polar (IR-inactive) mode. The situation 
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is very much different in the case of the IR-active mode. Indeed, the IR-active phonons induce electric fields which 
interact with the Josephson phase directly as indicated by the Josephson relation (Q). Such a coupling will lead to the 
dissipation of the Josephson plasma oscillations due to transferring its energy to the IR-active phonons. This effect 
has been studied quite well |lq| , and it may serve as a reference for assessing a significance of the IR-absorption in the 
case of the non-polar phonons. In order to make a proper comparison, we introduce some polar mode. It is a vector 
displacement field £'(x, t) in the barrier, which corresponds to the electric polarization density ea -3 £' (taken per 
unit cell of the volume a 3 ). Then, the interaction energy between the a.c. voltage V and £'(x, t) is simply the 
polarization energy 

H' m t=-^ J d 3 xE z (x,y,t)? z (x,y,z,t), (14) 

where E z is the electric field in the barrier, which is taken to have only one component along the normal to the 
barrier and is assumed to be uniform along this direction. Then, one finds E z = V/db , and E z can be related 
to the phase oscillations by means of the Josephson relation (^). The integration is carried out inside the insulating 
barrier of the width db- Thus, assuming that the barrier is thin, so that £' z can be taken uniform along the z-axis 
inside the barrier, we find eq.(|l4]) as 

H ' mt = ~2^? I tedyipfayrftzfaytt). (15) 

Thus, the polar mode interacts directly with the phase oscillations regardless of the presence of the stationary currents. 
Later we will employ the form ( [l5| ) as a reference for assessing the strength of the IR-absorption arising in the case 
of the non-polar mode described by eq. ( |l3|) . 

We consider two types of processes initiated by the term ( |l3| ) : i) oscillations £ can excite the J J plasma waves V 
at the same frequency, in accordance with the spatial profile of Jo(<£>o( x )); h) the incoming electromagnetic radiation 
excites the phase oscillations which convert into the oscillations £ directly at the same frequency. This constitutes the 
IR-absorption by the originally IR- inactive mode. The intensity of such an absorption is determined by Jo(^o( x ))- 
Accordingly, the stationary current distribution Jo(<yJo( x )) can be inferred. 



III. DIRECT CONVERSION OF NON-POLAR PHONONS INTO THE JOSEPHSON PLASMA 

OSCILLATIONS 

In this section, we consider the process i). To be specific, we choose y, z, t) as the plane wave A u exp[i(wi — 
kx)] + c.c., characterized by the frequency ui and by the wave vector k = (fc x , k y , k z ) as well as by some constant 
amplitude A^. The junction plane is taken at z — 0, and the insulating layer (the barrier) is taken as being of zero 
thickness. This approximation can be justified, if the actual barrier thickness db (which is about few nm) is much 
smaller than l/fc 2 , so that £ does not vary significantly inside the barrier along the normal. Hence, in what follows 
we will employ the form 

£ = A u exp[i(u>t — k x x — k y y)] + c.c. (16) 
at the location of the barrier in eqs.(|^),([T3|). Performing, then, the linearization of eq.(||), we obtain 

(-n 2 +^ 7 (x))^(x) - c 2 V 2 ^(x) +bvJ (<p (x))A,e- iik * x+k » v) = 0, (17) 

where we have introduced the notation ^(x) for the time independent amplitude of ^/>(x, i) = ?/> w (x)e' wt + c.c; 
f2 2 = uj 2 — ijuj; and 

2 dJo(yo) , 1Q s 

cVo 

is the square of the local Josephson plasma frequency. 

Eq.(p"7|) describes the electromagnetic linear response of the Josephson phase on the phonon field ( |l6| ) |19|. It is 
important to realize that the resulting a.c. Josephson voltage (|J) turns out to be given by the equilibrium Josephson 
currents Jo{tpa). This feature can be employed for the imaging of such currents. It should, however, be noted that 
the non- uniformity of Wj(x) ( |l8|) complicates solving eq. (fl7l) . In order to avoid this problem, we will consider the 
case 
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= |o/ >wj . (19) 

Taking such a limit, on one hand, allows omitting the term ~ uPj from eq.([l7|), which simplifies the solution 
significantly. On the other hand, eq. (|l9|) implies that the magnitude of the produced voltage is suppressed. A solution 
of eq.(0) for the uniform a.c. voltage V(t) = exp(iwi) + c.c., where V(t) — h/(2eS) J cPxip and the integration 
is performed over the whole area S — LW of the junction, with L and W standing for the length (along x) and 
the width (along y) of the junction, respectively, is 

V u = %MJ 0k , J ok - ( d 2 xJ (Mx))e- l{k * x+kyV \ Z £ (oj) = , (20) 



where Jok stands for the 2D Fourier component of Jo(tpo) along the GB; an explicit expression for v (see eq. 
has been employed. The solution (^0|) shows that the uniform a.c. voltage induced by the non-polar phonons incident 
on the JJ is proportional to the spatial Fourier component of the static Josephson current Jo(<^o( x )) at the wave 
vector (k x , k y ) (along the junction plane) of the incident phonon field. 

Let us estimate magnitude of the voltage ( p(j| ) in the non-resonant regime, that is, when the condition ( fl9] ) is 
satisfied. A typical extension of the current structure in the JJ is given by the Josephson penetration length Lj = 
y/hc 2 / (%iredJ c ) w 1/im for the critical currents J c w 10 5 A/cm 2 and d w 3 ■ 10~ 5 cm. Thus, the projection of the 
phonon wave vector along the junction plane (the k x ,y components) should be fen « Lj 1 w 10 4 cm _1 in order to 
produce the most efficient imaging. In this case, Jok ~ J c LjW for W <C Lj and L ^> Lj. Then, we choose 
bA u = 10~ 2 which corres pon ds to the maximum relative displacement (deformation) about 1CP 4 . We also employ 
the values c/c ~ 10 -2 yielding uij = c/Lj ss 10 12 s _1 . In order to satisfy the condition (|l9|), we choose 

il 2 w lOw 2 , which corresponds to 



|0| = y/\uj 2 - i 7 w| W 3^j. (21) 



When dissipation is negligible (large shunt resistance, so that 7 <C wj in eq.(|6|)), the condition (|2lj) can be satisfied 
by w> 10 13 Hz. While being in the range of the optical frequencies, so high u> is hardly achievable for sound. If, 
however, J c is in the range 1 — 10A/cm 2 , the plasma frequency becomes loj ~ 10 9 — 10 10 s _1 , and the condition 
( pl| ) can easily be fulfilled for sound. It should also be noted that a less demanding requirement for the sound can 
be imposed for the over-damped JJ, where 7 3> wj in eq.(^). In this case, already w 2 /7 < u> <C loj fulfills the 
condition (|2l]). This allows one to choose uj which is much less than the above values. For example, for 7 = 10 14 s~ 



and uij = 10 12 s 1 , the condition (21) is satisfied by u> = 100GHz, which is on the upper limit for the hypersound 
techniques pp]| . Finally, setting the length of the junction L m lOLj, we obtain the value w 10 _8 V in eq.([20|) 

Let us also justify the validity of the approximation neglecting eq.(S). As discussed above, taking eq.(f|) into account 
may result in significant deviations of the CPR from the standard one. This, however, does not occur if the condition 
( pd| ) holds. Employing the above values in eq.([i~ij) and choosing ujq ~ lo = 10 12 s _1 and v g — 10 5 cm/s, p « 6g/cm 3 , 
we find g « 3 • 10~ 3 <1 as mentioned above. 



IV. IR-ABSORPTION INDUCED BY THE STATIC JOSEPHSON CURRENTS 

Now let us consider the case ii), that is, the possibility of the coherent phonon field generation as a result of applying 
a uniform a.c. voltage exp(iu)t) + c.c. across the barrier. This voltage corresponds to the phase oscillations 

1>(t) = ^' + c.c, (22) 

following from eq.(||). These phase oscillations serve as a source (in eq. ([?])) for the oscillations of £ in the barrier, 
which are subsequently emitted into the bulk. In order to describe this effect, we employ eqs .([?]), (p^), and retain only 
the part ~ ip m the last term of eq. (0) . Accordingly, eq. (^) becomes 

I + c 2 £ - vl V 2 £ + ( ^M e"' + c.c) S(z) = 0, (23) 

where eq.(]2^) has been employed. Eq.(p3|) describes the linear response of £ on the external a.c. voltage. This can 
be interpreted as though the displacement field £ has become polar, so that it exhibits the IR-activity inside the 
JJ barrier. As a result, an additional dissipation of the external radiation occurs. The total power P{to) carried 
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away by the phonons, which are generated in the barrier and, then, are emitted into the bulk, can be calculated by 
means of solving eq. (|23|) and substituting the result into the expression for the phonon energy flux density —pVg£V£. 
Averaging this flux over time and integrating over the JJ area, we find 

P(w)=R (w)|VL| , i? (a;) = / 77^-^= (24) 



wIpU 2 , 7 (27r) 2 / -2/- 2 2\ 2 ' 

IP 9 17 v ; y Vg {lu 2 - lu 2 ) - q 2 



where lu > luq (for v 2 > ); and the integration is performed over the range |q| < <J v g 2 (lu 2 — ljq). The above result 
indicates that the frequency dependence of the absorbed power is determined by the spatial profile of the stationary 
current distribution. Thus, measuring the IR-absorption can be employed for obtaining information about J (<^o)- 
The quantity carries an obvious meaning of the a.c. conductance due to the phonons. 

It is useful to compare the absorption intensity (p^), caused by the emission of the non-polar phonons, with the 
intensity in the case when the phonons are polar, and the coupling is represented by eq.(|l5|). Such a comparison 
helps evaluating a significance of the effect under consideration. This is especially important because the technique of 
the IR-absorption in the J J is well established Jl^] , and therefore can be employed for detecting the effect described 
above. Thus, we assume that the displacement fields in the both cases are the same ( £ = In the case (|l5|), 
the last term in eq.(^3|) should be replaced by eV/{po?)5{z). The corresponding emission intensity P'{lu) can be 
calculated similarly to how eq.(p3) has been obtained. Comparing it with we find 



(25) 



where we have estimated the integrals in eq.(p4]) as J d 2 q\ J d 2 xJ a exp(iqx)| 2 w LjW\J c \ 2 , and has selected the 
emission area, from which P'(u>) is collected, as LjW. Obviously, in the case represented by eq.(p4|), LjW is 
the effective area occupied by Jo(ipo). Taking lu — 10 12 s _1 , a = 0.4nm, b — 10 2 /a and J c = 10 5 A/cm 2 , we find 
the ratio ( p5| ) as ~ 10~ 2 . For lu = 10 11 s _1 this ratio becomes close to 1. Thus, the IR absorption per unit area 
occupied by the equilibrium Josephson currents can be comparable to the absorption due to the IR-active mode. This 
conclusion makes detecting the effect described above by means of employing the technique E3] quite feasible. 

V. IMAGING OF THE STATIONARY JOSEPHSON CURRENTS BY SOUND 

The effect discussed in Sec. Ill gives a principal possibility of imaging the currents Jo(<^o( x )) by means of detecting 
uniform component of the a.c. voltage produced as a result of subjecting the barrier to a coherent time dependent 
field £. In this case, eq.(|o|) relates such a voltage to a single spatial Fourier component of Jq(lpq(x.)). The drawback 
of this is that, in order to keep such a relation valid, the frequency of the incoming phonon field should be high 
enough. As we will discuss below, this limitation can be circumvented by employing the reverse effect discussed in 
Sec. IV, that is, a coherent emission of the non-polar phonons induced by external electro-magnetic radiation. In this 



case, measuring the phonon field £ (in, e.g., far zone) which obeys eq.(23) would allow to achieve the same goal. We 
will discuss this in detail for sound, so that detecting a dynamical mechanical stress in the far zone could be employed 
for the imaging of Jo(</3o( x ))- 

In what follows the role of £ is played by the strain tensor Uij. We employ a simplest model of the dilatation 
interaction between « tJ and the Josephson currents, so that we take £ = Uu, and describe the dynamics in the bulk 
within the isotropic medium approximation. Thus, the energy explicitly dependent on ity is 

H d = J d? x [^ + \u\ + im% + if* «„], r, = ^p±S(z) , (26) 

where A, fi are the Lame coefficients; p stands for the total density; the source term ~ 77 is due to Hi nt (p"3|). 

Let us assume that small uniform phase oscillations ip(t) ( |22] ) are imposed externally by the incoming microwave 
radiation We represent the deformation field as u, = VjU, which determines Uij = (VjUj -I- \/jUi)/2, with 

u being some scalar field. We note that the transverse part of the deformations is decoupled in this simplified 



approach. Then, the equation for u , following from eq. (26), is 



^V 2 u = ^-V u B bA + cc, (27) 
ihiop 
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where v s — y/(A + 2p)/p stands for the speed of the longitudinal sound; eq.(^2|) has been employed; and no 
dissipation of sound is considered. Solution of eq.(27) can easily be found in terms of the Fourier harmonics u — 



M q e™"" lqx + c.c. as 

2eV u 



iTiujp(v 2 q 2 — uj 2 ) 



V« , (28) 



where ?7 q = J d 3 x exp(iqx)ry(x) denotes the 3D Fourier harmonic of 77 in eq.(p6|), and V w is taken uniform along 
the junction. As represented, eq.(p8|) is valid for an arbitrary configuration of the source r\. Thus, it is not limited 
by the specific form (2q) corresponding to the plane junction. 

Let us find the w-harmonic of the pressure field a = pv 2 uu = pvfV u far from the junction in the effective 
2D geometry. That is, we assume uniformity of all the quantities in the y direction corresponding to the thickness 
W of a slab containing the JJ plane z = 0. The other two directions x, z are unlimited. The location of the 
spontaneous currents determined by the solution of eq.(||) is restricted to some region of a typical size Lj (along x ) 
and W (along y). The distance r (along x, z directions) at which the stress is detected is considered to be much 
larger than Lj and W. Performing the inverse Fourier transform of eq.([2^) in the far zone, we find the w-harmonic 
of the pressure cr(x, z, t) = cr w (r, 9) exp(iujt) + c.c. as 

1 a\ e 9" /2 K, Vie~ zqr 
<W,#) = —7=7 ^Vq, (29) 

where the notation = uj/v s has been introduced, and q = g^n, with n standing for the unit vector in the direction 
of observation from the location of the currents, that is, r = (x, z) and r = rn, with n x — sin0, n z = cos 8 given 
by the angle of observation 9 with respect to the axis z. Eq.(]2^) can be employed for imaging of the spontaneous 
currents by means of detecting the mechanical pressure (stress) at the frequency of the external a.c. voltage far from 
the junction. 

For the purpose of estimating the magnitude of the stress, we assume that the equilibrium Joscphson current is 
uniform in the y-direction. Then, after employing the explicit expression for rj (|26|) in cq.(p9|), we find 



* u (r,6) = 7"_ y " v V / dxflf^MM*))- (30) 
2V27rw vr J 

This expression shows that, if the stress amplitude in the far zone is known for a sufficient range of the angles of 
observation 9 and frequencies uj, the stationary (spontaneous) Josephson current distribution can be restored by 
means of the spatial inverse Fourier transform. Later we will estimate the expected values of a. 

Now let us obtain an expression for the total electric power dissipated due to the emission of sound (see eq.(|2~i|) 
for optical phonons). As we will see, this power demonstrates a universal behavior ~ ui 3 as uj — > as long as the 
current distribution Jo is characterized by some finite vorticity. 

The dissipated power P(uj) equals to the total flux of energy 

P(uj) = -rW J dnaVii (31) 

carried away by the sound. Here the surface integral of the energy flux density is taken over the cylinder surface of 
the radius r ^> Lj and of the height W. A substitution of eqs.(E8h, (|30|) yields 



-iqr 



p{uj) = JrV)|K,| 2 , R-\uj) = ^ / <w 



dx e i(9 - sin9)x J (<Po( x )) 



2 



(32) 



It is important to note that the frequency dependence of i? _1 (w) turns out to be sensitive to the net vorticity. 
Indeed, let us consider the case when the J J contains some vortex characterized by the phase variation (fo(x) along 
the junction, so that the phase difference on x — ±00 is Sipo =/= 0. Then, taking into account eq.(^J) in eq.(p2[), we 
obtain 

R (uj)-R uj, R Q 2567r2e2v e d 2 p - (33) 

where the wavelength 2ir/q of the emitted sound is taken longer than the Josephson length Lj ( along which the 
variation of the phase typically takes place), so that the Fourier transform of Jq was taken as 
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Joa = — iq u sin 8 [ dxe~ 1 ^ sin e)x „ ff j g sm e6<p (34) 

v dx v 

in eq.(|3^), and we have employed the explicit expression for v in eq.(^|). Thus, the emission of sound induces the 
specific contribution R^iuS) ~ uj 3 (5ipo) 2 to the a.c. conductance of the J J in the limit ui — > 0. We note that in the 
absence of the vorticity ( 5ip = 0), eq.fl34|) yields Jn 9 ~ q 2 , and eq.(33) will change to become ~ w 5 . This 



result indicates that a presence of the spontaneous J J currents characterized by finite vorticity could be revealed in the 
frequency dependence (|33|), unless, of course, the low frequency quasiparticle effects dominate the JJ a.c. conductance. 

Let us estimate a possible magnitude of the effect (|33"|). For this purpose we choose the following parameters: Sfo = 
2tt; W=10- 3 cm; v s w 4 • 10 5 cm/s; p = 6g/cm 3 ; d = 3 • 10" 5 cm; b = 10 2 . This yields R~ x = 3- lO^Q-^w/GHz) 3 . 
It is important to note that the dependence (|33|) saturates at some value R m \ x at frequencies larger than a frequency 
t^max corresponding to the wavelength of sound comparable to the Josephson length Lj. For Lj = 1/xm, we find 
u-max ~ v s /Lj = 4-10 9 Hz. Accordingly, R^. rs 2T0" 3 O- 1 . This corresponds to ^^/(WLj) » 2-10 4 ft- 1 cm- 2 of 
the inverse effective shunt resistance per unit area of the J J. 

It is worth noting that measurements of i? _1 (w) do not allow restoring completely a spatial profile of the currents. 
As eq.(|32"|) indicates, i? _1 (w) is the quantity which is the average over the directions. Furthermore, the quasiparticle 
effects may dominate the a.c. conductance. In order to be able to restore Jo(y?o(x)), the mechanical stress amplitude 
( ^p| ) should be determined as a function of 8 and u>. Let us estimate its magnitude at distances r = 10 3 Lj far from 
the location of the currents and for q w Lj. The voltage amplitude in eq.(pO|) should be chosen in such a way that 
the linearized approach is valid. Thus, eV^ < Tllo. Then, the Fourier component of the current Jq q in eq.(pO|) can 
be evaluated as J c Lj, where a typical value of the critical current can be taken as J c w 10 5 A/cm 2 . Then, eq.(|30|) 
yields \a u \ w 0.1(w/GHz) 3 / 2 Pa for Lj = lpm. 

VI. ENHANCEMENT OF THE INTERCONVERSION BY EXTERNAL MAGNETIC FIELD 

It is worth noting that the discussed effects can exhibit a resonant-like behavior with respect to applied external 
magnetic field H y (along the J J width W ). Indeed, an external field greater than the lower critical field H c i w 
<& /{L,jd) [0, where <I>o = hc/2e is the magnetic flux quantum, leads to the periodic (along the junction length) 
structure of the Josephson currents due to the Josephson vortices. This periodic structure can contribute constructively 
to the induced voltage (^0|) and to the stress field (|30|), if the matching conditions are fulfilled. As a result, a significant 
enhancement of the corresponding values should be anticipated. Let discuss this. 

In the presence of H y > H c \ , the Josephson current will change its sign on the length Lh ~ ^o/{H y d) < Lj [ ^6[ . 
Thus, Lh describes the spatial period of the main harmonic of the current. Accordingly, as follows from eqs.(p0|), the 
induced voltage will exhibit the resonant increase when the wave vector fen (along the J J barrier ) of the incoming 
deformation field becomes equal to 2%/Lh, that is, for 

$ofe|| , , 

In this case, provided the vortex lattice is ideal, the Fourier component of the current becomes Jok ~ J C WL in 
eq.(]2(]). This implies the factor of L/Lj ^> 1 increase of the voltage magnitude if compared with the case of a single 
vortex considered in Sec. III. 

A similar situation occurs for the mechanical stress ([30j) generated by the incoming a.c. voltage. This stress will 
exhibit a resonant increase at certain angle of observation 8 given by the matching condition 2%/Lh = (uis'm8)/v s . 
Thus, for given frequency lu of the external radiation and the angle 8, the stress amplitude will resonantly increase 
at the value 

$^sin^ 



For long J J |16|, we find from eq.(|30|) that a u will increase in proportion to the total length L of the J J. This also 



constitutes the factors of L/Lj 3> 1 increase of the stress amplitude (30) in the given direction, if compared with 
the above estimates for the case of a single Josephson vortex. 

VII. DISCUSSION AND CONCLUSION 

The interconversion effect discussed above can, in principle, be expected to occur in traditional superconductors 
as well. Deformations should modulate the width db of the JJ barrier. Accordingly, the tunneling matrix clement 
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will change. In the superconducting state, this leads to the term (|T^), where the coefficient b is determined by the 
kinematics of the one electron tunneling. Assuming the exponential form ~ exp(— db/X) for the matrix element, 
where A stands for a typical overlap length of the electronic states on the both sides of the JJ, the estimate 
b f» db/X follows, if the dominant effect of the strain is the change db — * d&(l + u zz ). Obviously, in order to have 
b as large as chosen above b 10 2 , the critical current J c ~ exp(— db/X) should be practically zero. Thus, in 
the traditional JJ, b ~ 1 — 10, which lowers the above estimates of (|20| ) and (|3^) significantly. On the contrary, 
in the JJ made of the high-T c materials, the nature of the term ( |l3| ) is expected to be completely different. The 
proximity of the barrier to the superconducting state should result in anomalously strong sensitivity of the barrier 
to any external factor which shifts its equilibrium. In our estimates of b for the high-T c JJ, we have relied on the 
indirect evidence - the extremely strong sensitivity of the Josephson critical current to the misorientation angle in the 



GB J J. As suggested in Rcf. 11 1, the mechanical strain is the most plausible cause of this. 



As discussed in Introduction, the mechanical strain (stress) may have no significant direct effect on the barrier 
transport properties, if its role is reduced to only controlling the depletion of oxygen at the GB during the prepara- 
tion of the J J at high temperatures. If so, in the superconducting state, the external stress will produce no substantial 
modulations of the critical current, and the value of b will not be much different from the estimate for the traditional 
JJ. Thus, no the strong interconversion discussed in this paper is to be observed under such circumstances. Accord- 
ingly, either observing or not observing the discussed effect will help elucidate the role of the mechanical strain in 
suppressing the Josephson current. If the effect is observed, this will be a very strong indication favoring the concept 
of proximity of the barrier to the superconducting state fl2|,|l5| . 

If being strong enough, the dynamical effect proposed above, opens up a possibility for the high resolution imaging 
of the equilibrium Josephson currents in the high-T c materials. This is especially important for revealing the nature 
of the superconducting OP at the GB, and the mechanisms of the TRSB. The resolution is limited by the wavelength 
of the sound. Thus, for the frequencies higher than 10GHz, the resolution becomes better than 1/zm. 

In summary, we have suggested that in the Josephson junction, in which the critical current is sensitive to the 
mechanical (optical) deformations, a direct conversion of the non-polar phonons into the electromagnetic Josephson 
oscillations as well as the reverse process can occur, if this JJ contains spontaneous currents. We conjecture that 
these effects should be especially strong in the high T^-materials due to the proximity of the insulating layer of the JJ 
to the superconducting state. We propose that the effect of the interconversion can be employed as a crucial test for 
the role of deformations in the superconducting properties of the GB JJ made of the high-T c materials. This effect 
can be used for imaging of the Josephson current spatial distributions. 
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